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The Davey-Stewartson (DS) equations with a perturbation term are pre- 
sented by taking a fluid system as an example on an uneven bottom. Stability 
of dromions, solutions of the DS equations with localized structures, against 
the perturbation is investigated numerically. Dromions decay exponentially 
under an effect of the perturbation, while they travel stably after the effect 
disappears. The decay ratio of dromions is found to have relation to velocities 
of dromions. The important role played by the mean flow, which acts as an ex- 
ternal force to the system, is discussed. These results show that dromions are 
quite stable as a localized structure in two dimensions, and they are expected 
to observed in various physical systems such as fluid or plasma systems. 
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I. INTRODUCTION 



Nonlinear problems in one spatial dimension have been studied intensively and there 
are a lot of fruitful results from these researches. One of them is the study on nonlinear 
integrable equations, such as Korteweg - de Vries equation or nonlinear Schrodinger equation, 
and soliton is derived as a solution of those equations. Solitons are localized structures 
with exponentially decaying bases, stable under interactions between each other, and have 
properties both of wave and of particles. Many problems have been discussed for these 
solutions from the viewpoints of stability under many kinds of perturbations [|T]-|^| . 

On the other hand, in the systems in two spatial dimensions or more, few equations are 
known to have localized solutions. The Kadomtsev-Petviashvili (KP) equation is one of the 
equations in two dimensions which have soliton solutions, but these solitons localize only 
in one direction. None of the known solutions of the KP equation other than the soliton 
solutions, such as lump solution or the resonance solution |J, localize in two dimensions with 
exponentially decaying base either. The solutions of the KP equation cannot be recognized 
as localized structures in two dimensions. 

Apart from the usual soliton solutions, we can nominate dromion H?HT0|] , which is a exact 
solution recently found in the Davey-Stweartson 1 (DS1) equations |II]]I^fl : 



as a localized structure in two dimensions. For this solution, the envelope of \A\ has a 
localized structure in two dimensions with a base which decays exponentially. 

Dromions have quite different properties compared to solitons, although both of them 
are localized spatially. For example, if two dromions have interaction between each other, 
the number of pulses generally may not be conserved. Hence we cannot infer behaviors of 
dromions from the analogies of those of solitons. Moreover, detailed discussions on stability 
against perturbations of dromions are quite few compared to those on solitons in one dimen- 
sional space. For these reasons, numerical analyses on the behaviors of dromions against 
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various perturbations are quite necessary. 

In this paper, our purpose is to take perturbations to the DS1 equations into account 
and clarify how they affect on the dromion solutions. We consider perturbation terms due 
to inhomogeneity of the bottom of the fluid in two dimensions. This paper is organized 
as follows. In the next section, we derive the DS equations with a perturbation term and 
consider the effect on dromions. In §3, detailed studies on the stability in typical situations 
is presented by numerical analyses. Discussions and concluding remarks are presented in 
the final section. 



II. MODEL EQUATIONS 

In the previous paper JIB], we presented the perturbed DS1 equations by using the 



reductive perturbation method. First we briefly summarize the derivation in order to see 
the meaning of the perturbation term, and then discuss the effect of it. 

We consider a system of irrotational incompressible fluid on an uneven bottom in three 
dimensions. The velocity potential will be expressed as </>(x, y, z, t). Then the basic equations 
are the Laplace equation 

A0 = O, (-H<z<(), (2.1a) 
and the boundary conditions at the surface and the bottom of the fluid. These are 

H x <j> x + H y <t) v + <t> z = 0, (z = -H), (2.1b) 

<fiz - Ct - <fix(x - <Py(y = 0, (Z=C)> ( 2 - lc ) 

<f> t + g( + (V<i>) 2 /2 = T/Rp, (z = (), (2.1d) 

where we have expressed the density of fluid as p, the bottom as z = —H(x,y,t) and the 
surface asz = ((x,y,t). The quantity g is the acceleration due to the gravity, T the surface 
tension and R the curvature of the surface, which is given by 



Now we are going to derive nonlinear amplitude evolution equations for the velocity potential 
by means of the reductive perturbation method |[Tl] . [T^ , p^5|] . The quantity H is not a constant, 



but we consider that its variation is slow and small in space. We also assume the spatial 
variance of ( is sufficiently slow compared with the characteristic length of the localized 
structures. Then let us consider the deviation of (0, (, H) around the stationary values 
(0,0, H^), where H<® is a constant. We expand them as 

oo n 

0=E E ^ B) & V, *, r)e n e^-^\ (2.2a) 

n=l l=—n 



oo n 



C = E E Cf(^V,z,T)e n e^-^, (2.2b) 

n=l i=—n 

oo 

H = H {0) + E H {n) (£, T], r)e n . (2.2c) 

n=l 

We have introduced expand parameter e and new independent variables £ = e(x — vt), 
7] = ey, z = z and r = e 2 t, where v is the group velocity. We assume relations = (j)^* 
and (}™} = Q from the reality of the physical quantities. The asterisk indicates complex 
conjugate. Substituting these into and writing K = k 2 T/pg and a = tanh kH (°> , we 
can get the dispersion equation 

u 2 = gka(l + K), 

and relations between expansion amplitudes (j)^ and in ( |2.2a|) and ( 2.2b|) . The main 



results up to 0(e 2 ) are listed in the Table Q. In the next order, we have a set of equations 
whose dependent variables are ^g 1 "* and (f)[ . we introduce new dependent variables A and 
Q instead of (f)^: 

A = A U 

Q = v{A 0( [gH {0) - v 2 ](l + Kf/k 2 + \A x \ 2 [v{l + o 2 )(l + Kf + (1 + K) - 2Ku/k)}, 

where A and A 1 are listed in the Table |I[ Applying suitable scaling and 45° rotation of the 
coordinates, we have 

iA t + % + A m - 2\A^A + A(Q^ + Q v ) = -iJ(f , rj)A, (2.3a) 
Q^ = {\A\\ + m\, (2.3b) 



when the coefficients of the amplitudes satisfies a certain relation. The function J(£, rf) is real 
and proportional to , with positive proportional coefficient. These are the DS1 equations 
with an extra term in the right-hand side of ( |2.3a| ), which comes from the inhomogeneity of 
the system. 

As we have discussed in the Ref. [[TjJ, we take proportional to tanh(p^ + qr)) (p and 
q are real constants). This means there is a mild step into one direction on the bottom. In 
this simple calculation shows that the perturbation term becomes 

iJ(£, rj)A = iasech 2 (p£ + qrj). 

The quantity a indicates the degree of inhomogeneity, and p and q determines the width 
and the direction of the step. 

We consider a traveling wave solution, and set A ~ F(t)e i( - k ^+ k yv-^) . This wave travels 
away from the origin in the first quadrant asymptotically. Since we have assumed that H 
varies slowly in space, we can neglect higher order spatial derivatives of Ignoring the 

nonlinear term in ( 2.3a ), we have u = (k 2 , + k 2 ) and F = exp[— a J(£, rj)t] ]T3|] . For the 



example we have introduced, since oc tanh(jo£ + qr)), the perturbation term will be 
proportional to sech 2 (p£ + qrj). This shows that the perturbation term results exponential 
decay or amplification of the pulses monotonously. 

III. RESULTS OF NUMERICAL ANALYSES 



In this section, we shall first present the localized solution of (|1 . 1|) used as initial-boundary 
conditions, and then show the results of the numerical analyses of the perturbed DS1 equa- 



tion ( |2.3j ). As we have pointed out in the Ref. |17]], it is inevitable to analyze the DS1 equation 



numerically in the case of dromion solutions, because of inexistence of the Hamiltonian and 
complexity of them. 
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A. Localized structures as an initial-boundary-condition and the numerical method 



Since we are going to focus on the stability of a dromion against perturbations, hereafter 
we restrict ourselves to the initial condition in which only one localized pulse appears. First, 
we shall show the explicit form of the one dromion solution for the DS equations ( |1 . 1| ) . 
Setting the unknown variables as [ITH] 



A — G/F, U = 2(logF) xx , V = 2(logF) yy , 



(3.1a) 



where the functions F and G are 

F = 1 + exp(77i + 77*) + exp(77 2 + Vb) + 7 ex P(^i + Vi + V2 + 

G = pexp(?7i + T] 2 ). 
Writing x — x\ and y = x 2 , we have defined the quantities rji and r\ 2 as 



+ ik 



U) 



Xj + 



a® + id? 



t, (j = l,2). 



From (|1.1| ) and (|3.1|) , we have the following relations between parameters 
\p\ = 2 V /2^ 1) ^ 2) ( 7 -1), 



Ej=l,2 



ki^ 2 - k 



0)2 



The boundary conditions for U and V are given by 



U 



8k^ 2 exp(r/i + r/*) 
y =-co [1 + exp(rji + T}*)] 2 



V 



8kP 2 exp(r/ 2 + T}%) 
x=-oo [l + exp(r/ 2 + r/!)] 2 



(3.1b) 



(3.1c) 



(3.2a) 



(3.2b) 



From ( Hm-3 ) and ( ^Tc|) , we can find that there are five free parameters in this solution. The 
quantity 7 determines the amplitude of the single pulse of \A\; The parameters k^ is the 
width of the pulse of \A\ in the x direction (when j = 1) and y direction (j = 2), while k\^ 
give the velocities in each of the directions. Hereafter in this section, we are going to use 
the initial and boundary conditions given by these ( [3.1|) and Q3.2|). 
Next we rewrite the equation to be analyzed as: 

iAt + A xx + A yy + (U + V)A - iH(x, y)A = 0, 



U y = {\A\ 2 ) X , V x = (\A\ 2 ) y , 



(3.3) 
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which is derived from (|2.3| ) by introducing the new variables U and V by 

U = Q X -\A\\ V = Q y -\A\ 2 . 

As we have considered in the previous section, the perturbation term is chosen as 
—iaAsech 2 (n x x + jj y y), where a is the degree of perturbation, and fx x and \i y determines 
the width and the direction of the step on the bottom. 

Now we have to choose the five parameters in the one dromion solution, 7, and k^ 
(j = 1,2), and the three which determines the perturbation term, a, \i x and \i y . Since we 
are going to consider a perturbation for an identical dromion, we use fixed values for 7 and 
kr \ (j = 1, 2) all through the simulations. The values are 7 = 3, jfeC?) = 0.6, (j = 1,2). We 
choose suitable values of a, fi x and fi y for each case of the simulations. 

Finally, we show the numerical method briefly. The space derivatives in (|3.3| ) is performed 
by using the psudospectral method |L9| with periodic boundary condition. Integration of 



the first equation of ( |3.3| ) in terms of t is performed by the Burilsh and Store method 
Tj| and the fourth order Runge-Kutta method with appropriate accuracy of adaptive step 



size control. The second equation of ( |3.3| ) is calculated by the fourth order Runge-Kutta 
method with boundary conditions which are presented in ( |3.2b[ ). When the midpoint value 
between meshes is necessary, we use cubic spline. The detailed discussion of the scheme of 
the calculation is described in the ref. [TtJ . 



B. The case when the dromion traverse the step 

The first analysis is for the case of fc| = kf^ = 0.5, fi x = —jj, y = 0.4. In this case, the 
contours along the step on the bottom runs parallel to the line x = y, which is chosen as the 
trail of the single dromion in the initial state travels along. This means that the dromion 
undergo the effect of perturbation constantly. The discussion in the previous section tells 
us that the height of the peak of the variable \A\ is considered to decay exponentially. We 
studied the time development of the height of the peak of \A\ for the various values of a, 



S 



a = 0.01, 0.05, 0.1, and 0.4. The results of this simulation is presented in Fig. |l[ The data 
of the case a = 0.4 deviates from a straight line with the lapse of time. This is because the 
pulse decays sufficiently (less than 1% in height), and hence it cannot be distinguished from 
the ripples distributed other places on the plane. We can immediately see that the result is 
consistent with the discussion above, and can also corroborate the accuracy of the scheme. 

C. The case when dromion goes across the step 

In this subsection, we set Kp = kf^ = 0.5, [i x = /i y = 0.4. This means that the initial 
dromion travels across the step. The dromion feels the effect of perturbation only when it 
passes near the origin. If dromions are not stable structure, they are supposed to break 
into small fractions due to the effect of perturbations, even after the effect vanishes. On the 
other hand, if they are stable, the opposite will occur; A localized pulse travels stably after 
disappearance of the effect of perturbation, though the height of its peak changes by the 
effect. 

The numerical analyses of the time development of \A\ is performed both for positive 
and negative values of a. In Fig. ||] and |3|, the time developments of the height of the peak 
is shown. The height of the dromion changes near the origin, only where the perturbation 
effect appears, but it keeps the same value after the dromion passed there. We displayed 
in Fig. f| the surface pictures of the main flow before and after the pulse get through the 
perturbation area. The dromion keeps a feature as an localized pulse, though it experienced 
perturbations. 

Next we select various values of the velocity of dromion and investigate how the effect 
of perturbation affects on the pulse. We consider a relation between velocities of dromion 
and decay ratio, the height of the peak of \A\ after dromion go through the perturbation 
normalized by its initial value. We choose a fixed value of a = 0.1, and different values of 
velocities = k^ = N/16, N = 8, . . . , 12. The relations between logarithm of the decay 
ratio and the inverse of the parameter fe t - , (j = 1, 2) is presented in Fig. |5|. This shows that 
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there is a linear relation between them. 

This result is explained as follows. As we have discussed in the previous section, a wave 
with small amplitude decays as F ~ exp(— a J(x, y)t). In this case, F(t) is the maximum 
value of \A\ and J(x,y) = sech 2 [/i(x + y)\. Considering the dromion traveling along the 
line y = x with velocity v, we have dF\ x ^ x+( [ x = —aFv~ 1 sech 2 (2fix) dx. Integration of this 
equation gives the function F to be F ~ exp[— a{2^iv ) _1 / <i£sech 2 £]. Then the decay ratio, 
D = F(oo)/F(— oo) satisfies a relation InD ~ — (a/yu)f _1 , and this shows a linear relation 
between InD and l/t>. 



IV. CONCLUDING REMARKS 

We have investigated time developments of surface waves described by the DS1 equations 
with perturbation numerically. The main results of the calculations tell that dromions decay 
or are amplified when they feel perturbation, and they travel stably after the effects of 
perturbations disappear. In addition, if the time during which the dromions undergo the 
effects of perturbation becomes longer, the degree of decay will be larger. Taking these 
numerical results into account, we can conclude that dromions are quite stable structures in 
two-dimensional space. 

Now we would like to consider the process under which the dromions travel stably in 
spite of the existence of perturbations. As we have discussed in the previous paper, the 
mean flow U and V play the important roles in ( |1 . 1\) or ( |3.3| ). They behave as attractive 



potentials of the Schrodinger equation, and hence preserve the localized structures of the 
main flow. Because the cross point of the peaks of U and V is the most attractive point 



in the system ||17|| , and it gathers structures of the main flow. The stability of the dromion 
solutions owes much to the mean flow. 

We can say that, as long as attractive points generated from the mean flow U and V 
exists, localized structures of the main flow will keep on existing. The stability revealed in 
this paper will assure observations of dromions in existing system, such as fluid dynamics 
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and plasma systems, and this will be a future problem. 
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FIG. 1. A semi- log graph of time development of the highest part of \A\ of one dromion for the 
perturbed DS equations when the dromion traverse the step. The height is normalized by that of 
the solution with no perturbation. The values of a are 0.01, 0.05, 0.1 and 0.4. The peak decays 
exponentially. The reason of the deviation from a straight line in the case of a = 0.4 is that the 
pulse gets as low as ripples distributed other places on the plane. 

FIG. 2. The developments of the highest part of \A\, when the dromion goes across the step. 
The values of a are 0.01, 0.05, 0.1 and 0.5. Since a is negative, the dromion goes into the direction 
where the depth of the fluid grows deeper, and undergo decay. The modification in height occurs 
when the dromion passes near the origin (t = 0), where the perturbation term has finite value. 

FIG. 3. Time developments of highest part of \A\. The motion of the dromion is the same as 
those presented in Fig. ||. The values of a are —0.01, —0.03 and —0.05. and this means that the 
dromion climbs up the step. The dromion is amplified, and the change of height occurs around the 
region where the perturbation exists as in Fig. ^. 

FIG. 4. Typical surface figures of the envelopes of the main flow, when the dromion goes across 
the step and decays. The parameter is a = 0.2, and the time is (a) t = —5.5, (b) t = —1.0 and (c) 
t = 4.5. Though the height of the peak is lowered, the pulse keeps the characteristic as a localized 
structure after the effect of perturbation. 

FIG. 5. The relation between velocity of dromion and decay ratio. The vertical axis is the 
logarithm of the decay ratio and the horizontal one inverse of k± , (j = 1,2). The parameter is 
a = 0.1. The straight line is drawn by the method of least square. The dromion which travels 
faster undergo less decay. 
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TABLE I. Main results of the expansion amplitudes in (|2.2[) from the reductive perturbation 
method. The variables Aj (J = 0,1) and Bj (j = 0,1,2) are functions of £, r] and r, and Z is 
defined by Z = z + H^ ' . The quantity v is the group velocity du jdk and c is the phase velocity 



ijjjk. The values with negative t can be derived from the relations 



An) 



and C 



(n) 



(n)* 



d n) (e,r/,r) 



(1,0) 
(1,1) 
(2,0) 
(2,1) 
(2,2) 



^i(£> ??> r ) cosh(fcZ) 
cosh(fe^ (0) ) 

Bo(e,»?,r) 

i?i (£, ^, r ) cosh(/cZ) — iA\^Z sinh(A;Z) 

cosh(/ci? (0) ) 
B 2 ^,v,r) cosh(2kZ) 



cosh(2kH 



o 

iuiA\ 
9(1 + K) 

vA oi 
~9~ 



2cK 



g(l + K) 



+ K 



k 2 A\ 



2koB\ 

25(1 + ^) ' ' w (l + <r* 



+ i- 
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